Magnetic order and spin-waves in the quasi-lD S=l/2 antiferromagnet BaCu2Si207 
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Elastic and inelastic neutron scattering were used to study the ordered phase of the quasi-one- 
dimensional spin- 1/2 antiferromagnet BaCu2Si2 07. The previously proposed model for the low- 
temperature magnetic structure was confirmed. Spin wave dispersion along several reciprocal-space 
directions was measured and inter-chain, as well as in-chain exchange constants were determined. 
A small gap in the spin wave spectrum was observed and attributed to magnetic anisotropy efi'ects. 
The results are discussed in comparison with existing theories for weakly-coupled quantum spin 
chain antiferromagnets. 



PACS numbers; 75.25.-f-z, 75.10.Jm, 75.40. Gb 
I. INTRODUCTION 

The investigation of low-dimensional magnets has been 
a very active research field for the last two decades. First 
conceived as simple spin models, low-dimensional mag- 
nets revealed a far richer physical behaviour than their 
more conventional three-dimensional (3D) counterparts, 
due to the importance of quantum fluctuations. Particu- 
larly interesting is the one-dimensional (ID) Heisenberg 
antiferromagnetic model (ID HAF), for which the Hamil- 
tonian is written as 

chain 

H = jY,S,- , (1) 

i 

J being the exchange coupling constant. The ground 
state of the ID S=l/2 HAF is a spin-singlet and can 
be characterized as "marginal spin liquid" : while long- 
range order is absent, spatial spin correlations decay ac- 
cording to a-|power-law, and thus represent quasi-long- 
range order.EJ The excitation spectrum is isotropic and 
gapless, and is described in terms oiyS =1/2 elemen- 
tary excitations known as "spinons".0 Physical excita- 
tions are composed of pairs of spinous, which gives rise 
to a 2-spinon-continuum.d This behavior is in stark con- 
trast with that of 3D spin systems, that have long-range 
order in the ground state ("spin solid"), and where the 
spectrum is dominated by single-particle 5 = 1 excita- 
tions (spin waves). The pure one-dimensional HAF is 
of course a physical abstraction since spin chains in real 
materials are always at least weakly coupled. To make 
a virtue of necessity, weakly coupled antiferromagnetic 
chains offer the opportunity to study the cross-over from 
the quantum spin dynamics in a single chain to semi- 
classical spin waves in 3D systems. 



BaCu2Si2 0p| was recently recognized as an almost ideal 
model systcniQ for experimental investigation of this di- 
mensional cross-over. The silicate crystallizes in an or- 
thorhombic crystal structure, space group Pnma, and the 
lattice constants are a = 8.862(2) A,b= 13.178(1) A and 
c = 6.897(1) A.Q The Cu^^-spins are subject to a strong 
super-exchange interaction via the 0^~-ions, producing 
weakly-coupled antiferromagnetic chains that run along 
the crystallographic c-axis.Q Correspondingly, the mag- 
netic susceptibility shows a broad maximum at around 
150 Ka and its temperature-dependence follows the theo- 
retical Bonner-Fisher curvdj with an intrachain coupling 
J = 24.1 meV. The nearest- neighbor Cu-Cu distance 
within the chains is equal to c/2 (for a schematic view of 
the crystal structure see Fig. 1 in Tsukada et aln) . Weak 
interactions between the chains result in long-range an- 
tiferromagnetic ordering at Tjv = 9.2 K, as obser\{ed by 
bulk susceptibility and specific heat measurements. □ Pre- 
liminary neutron diffraction studies provided an estimate 
for the ordered magnetic moment at low temperatures: 
Too = 0.16/LtB.a In the ordered state the spins are aligned 
along the crystallographic c-axis (chain axis). Early in- 
elastic work suggested that the ratio of inter-chain to in- 
chain interactions J'/J in BaCu2Si207 is intermediate 
compared to that in such well-known = 1/2 systems 
as KCuFa (a more 3D-like material) and Sr2Cu03 (an 
almost perfect ID compound), as shown in Fig. |l]. 

In a recent short paper we reported studies of the inter- 
play betweerusingle-particle and continuum dynamics in 
BaCu2Si207.a For a better quantitative understanding of 
this behavior a detailed knowledge of both in-chain and 
inter-chain interactions in this system is required. The 
present article is a detailed report on neutron scattering 
characterization of this material. We focus on the 3D as- 
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pect, namely long-range magnetic order and spin wave- 
like excitations. From measurements of the spin wave 
dispersion along different reciprocal-space directions we 
determine all the relevant parameters of the system, such 
as the strengths and geometry of inter-chain coupling and 
magnetic anisotropy. 
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FIG. 1. The ordered moment and Tm / J for the known 
weakly-coupled antiferromagnetic Heisenberg spin-1/2 chaki 
compoundsi as a function For Sr2Cu03tj, 

Ca2Cu03tZl and BaCu2Ge2 07l3til, |J'j was estipated using 
the predictions of the chain- mean field model,li3 "'^^jf^' fac- 
tual measured values were used for Tjv. For KCuFsErta \J'\ 
was deduced from spin wave bandwidth perpendicular to the 
chain ^is. The lines are predictions of the chain mean-field 
modeiia. 



II. EXPERIMENTAL 

Single crystals were grown using the floating-zone tech- 
nique from a sintered polycrystalline rod. The crystals 
were carefully checked for twinning because the lattice 
constants along the a- and c-axis are very similar. The 
diffraction measurements were performed with a small 
sample with dimensions 3 mm^ to exclude sizable extinc- 
tion effects. The samples used in the inelastic scattering 
experiments were cylindrical in shape, 5 mm in diameter 
and 50 mm in height. Up to three single crystals were 
co-aligned for the different experiments. 

The neutron scattering experiments were performed 
on five different instruments. The two-axis D23 diffrac- 
tometer at the Institut Laue-Langevin (ILL), Grenoble, 
France, was used to investigate the ordered magnetic 
structure at low temperatures (Setup I). D23 is a double- 
monochromator thermal-neutron diffractometer with a 
lifting detector arm. The experiment was performed 
with graphite monochromators and an incident energy 
Ei = 14.7 meV. 

Spin wave dispersion perpendicular to the chain axis 
was measured using two cold-neutron triple-axis spec- 



trometers: the TASP spectrometer at Paul Scherrer 
Institut (PSI), Villigen, Switzerland (setup II), and 
the IN14 spectrometer at ILL (Setup III). The zone- 
boundary energy was measured using the thermal neu- 
tron triple-axis spectrometer IN22 at ILL (Setup IV). 
Preliminary measurements were also performed at the 
National Institute of Standards and Technology Center 
for Neutron Research (NCNR). 

The experiment on TASP was performed with fi- 
nal neutron energy fixed at = 8 meV. PG(002) 
reflections were used for monochromator and ana- 
lyzer. The sample was mounted with its (/i, h, I) and 
{h, 0, 1) crystallographic plane in the horizontal scattering 
plane. The supermirror neutron guides gave a source-to- 
monochromator coUimation of 70' both horizontally and 
vertically. The measurements in the {h, 0, Z)-plane were 
performed using 80'-80' collimators before and after the 
sample and distance collimation after the analyzer. For 
the measurements in the {h, h, Z)-plane, the use of 2 cm 
wide slits cut the beam divergence before and after the 
sample down to effective values of 47'-60'. The energy 
resolution was 0.38 meV and 0.6 meV for the measure- 
ments in the {h,h,l) and (/i, 0,Z) crystallographic plane, 
respectively, as determined from the full width at half 
maximum (FWHM) of the quasi-elastic peak. 

For the measurements using IN14 and IN22, PG(002) 
reflections were used for monochromator and analyzer. 
The sample was mounted with its (0, k, I) crystallo- 
graphic plane in the horizontal scattering plane. Final 
neutron energy was fixed at i?/ = 3 meV and Ef = 
30.5 meV, for the IN14 and IN22, respectively, with neu- 
tron guides defining the source-to-monochromator coUi- 
mations of 30' and 60', respectively. 40'-40' collimators 
were used in both setups before and after the sample, and 
the use of slits cut the effective collimation of the IN22 
setup down to 27'-40'. The energy resolution (FWHM) 
was 0.26 meV and 1 meV for the IN14 and IN22 measure- 
ments, respectively. Be- or pyrolytic graphite filters were 
used after the sample to eliminate higher-order beam con- 
tamination. 

The time-of-flight chopper spectrometer MARI at the 
ISIS Facility, Rutherford Appleton Laboratory, Oxon, 
United Kingdom, was used for the measurement of the 
magnetic excitations in a wide range of wave-vector 
transfers along the chain (Setup V). The measurements 
were performed using a chopper-monochromated incident 
energy Ei — 30 meV allowing the measurement of the ex- 
citation spectrum up to 25 meV. The energy resolution 
was about 0.35 meV, as determined from the FWHM of 
the quasi-elastic peak. The resolution of the wave-vector 
transfer was given by the size of the sample, the detector 
and their respective distance, and it was typically 0.05 
r.l.u. along the chain-axis. 
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III. MAGNETIC STRUCTURE 

To determine the spin structure in the ordered phase, 
48 magnetic Bragg intensities were measured at T = 
1.5 K using Setup I. The main comphcation in the ex- 
periment was that all magnetic reflections appear on 
top of strong nuclear peaks. The following measure- 
ment procedure was implemented. For each reflection 
the nuclear intensity was measured in a rocking curve at 
T = 10 K > Tm- The peak intensity was then measured 
with high statistics both above and below Tjv, and the 
diff'erence was attributed to magnetic scattering. The in- 
tegrated magnetic intensity was estimated from the mea- 
sured intensity ratio and the integrated nuclear intensity. 
The magnetic structure was found to be totally consis- 
tent with that suggested in Tsukada et al.Ja The ordered 
moment at low temperature mg = 0.15(l)/LtB is parallel 
to the crystallographic c-axis. Relative nearest-neighbor 
spin alignment is ferromagnetic along the a-axis and an- 
tiferromagnetic along the 6- and c-axes, respectively. 
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FIG. 2. Temperature dependence of the (0 1 l)-reflection 
between 1.5 and 15 K. The scattering intensity above 
Tn = 9.02 K comes from the allowed nuclear Bragg reflection 
at (0 1 1). Below Tn, the intensity increases with decreas- 
ing temperature. The solid line is a fit to the data with the 
critical exponent /3 = 0.25 as explained in the text. 



The temperature dependence of a few magnetic reflec- 
tions was measured in the range T = 1.5-15 K as shown 
in Fig. 1^ for the (0 1 l)-peak. Below the ordering temper- 
ature T/v, the Bragg intensity increases with decreasing 
T. The measured peak intensity / is a sum of nuclear and 
magnetic contributions, the latter being proportional to 
the square of the sublattice magnetization m. The data 
were fit to a power-law form: 



I{T) = lo 



( Tn-T 



N 



2/3 



TV for T < Ti 



N ■ 



(2) 



where /q is the magnetic intensity at T = K, /3 is 
the order-parameter critical exponent, and N is the (nu- 



clear) scattering intensity assumed to be T-independent. 
A good fit is obtained with fi = 0.25(5) and T/v = 
9.0(0.05) K. The value of (3 is clearly below the critical 
exponent for a 3D antiferromagnet, and is in good agree- 
ment with the critical indexes nieasured in the quasi- ID 
materials KCufJ and SraCuOaN. 



IV. DISPERSION OF MAGNETIC EXCITATIONS 

A. Dispersion along [/i, 0, 1] and [h,h,l]. 

Spin wave dispersion at the ID AF zone-center q\\ = tt 
was measured along the [1,0,0] and [1,1,0] reciprocal- 
space directions at T = 1.5 K using Setup II. Typi- 
cal background-subtracted constant-Q scans are shown 
in Figs. ^ and 0. The background was measured at 
Q = (/i, 0,1.25) (/i=0-l) and Q = (/i,/i, 1.25) (/i=0-l), 
where no magnetic scattering is expected in the relevant 
energy transfer range, due to the steep dispersion along 
the chain axis. All measured scans show typical spin 
wave peaks with a sharp onset on the low-energy side, 
and an extended high-energy resolution tail. 
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FIG. 3. Constant-Q scans measured at T = 1.5 K for three 
different wave- vector transfers on the [/i, 0, 1] reciprocal-space 
rod (Setup II). The background has been subtracted as de- 
scribed in the text. The solid lines represent global fits to the 
data using a model SMA cross section function. 



The data were analyzed using a model cross-section 
written in the single- mode approximation (SMA), as will 
be discussed in detail in section The cross sec- 

tion was numerically convoluted with the spectrometer 
resolution function, calculated in the Cooper-Nathans 
approximation.E3 Global fits to all data, as well as fits 
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to individual scans were performed in the energy trans- 
fer range 0-4.5 meV. As shown in Figs. |^ and |[ the 
model used gives an accurate description of the experi- 
mentally observed peak shapes. The transverse spin wave 
dispersion relation deduced from the global fit is shown 
in Fig. ^ (solid lines), where symbols represent excitation 
energies obtained in fits to individual scans. 
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FIG. 4. Constant-Q scans collected at T = 1.5 K for three 
different wave-vector transfers on the [h, h, 1] reciprocal-space 
rod (Setup II). Background and lines as in Fig. 



B. Dispersion along [0, k, 1] 

Dispersion along the [0, fc, l]-direction was measured 
using the high-resolution cold-neutron Setup III dX T — 
1.5 K. Constant-Q scans were collected at wave vectors 
with k in the range Q < fc < 3. Some of these data were 
presented elsewhere^ Scans at the 3D AF zone-centers 
(0, 1, 1) and (0, 3, 1) (Fig. |) reveal the presence of a small 
gap in the spin wave spectrum. In fact, there appear to be 
two separate modes with slightly different gap energies. 
These two branches correspond to different polarizations. 
The relative intensity of the lower-energy excitation at 
Q = (0, 3, 1) is clearly smaller than at Q = (0, 1, 1). Due 
to the intrinsic polarization dependence of the magnetic 
neutron cross section, the observed behavior is consis- 
tent with the higher-energy mode being polarized along 
the a-axis of the crystal, and the lower-energy excitation 
polarized along h. This consideration is made quantita- 
tive by a global fit to all measured scans using an SMA 
cross section with built-in neutron polarization factors 
for the two branches (solid lines in Fig. l5|) . 
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FIG. 5. Constant-Q scans collected at T = 1.5 K for 
Q = (0,0,1), Q = (0,1,1) and Q = (0,3,1) (Setup III). 
The solid line is a global fit to the data as described in the 
text. Shaded areas indicate partial contributions of two spin 
wave branches with different polarization. 



The spin wave dispersion relation along the [0, l\ 
reciprocal-space direction is shown in Fig. ^ (solid lines). 
Here symbols represent fits to individual scans. 
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FIG. 6. Dispersion of the spin-wave excitations at 
T — 1.5 K along three different crystallographic directions 
perpendicular to the chain axis. The open and closed circles 
show the excitation energies extracted from fits to individual 
scans. The open circle indicates that the measurements were 
actually performed at (0, 1 — fc, 1), ( — ft, 0, 1) and {—h, —h, 1), 
respectively. Error bars are smaller than the symbol size. The 
solid line is the dispersion determined in a global fit to the 
data using Eq. ^. 
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C. Dispersion along the chain axis 

The dispersion of magnetic excitations along the chains 
was previously determinedrin constant-energy scans up to 
16.5 meV energy transfer.u As part of the present work 
we performed additional measurements of this type us- 
ing the time-of-flight technique (Setup V). The incident 
neutron energy at the MARI spectrometer was set to 
30 meV and the chain axis was aligned nearly perpen- 
dicular to the incoming beam direction. The scattered 
neutrons were counted by a large array of detectors, ar- 
ranged in a half circle with radius 4 m vertically below 
the sample. The energy transfer was determined from the 
detection time. This allowed a simultaneous measure- 
ment of the neutron scattering cross section on a two- 
dimensional surface in (Q, a;)-space. The projection of 
the measurement surface onto the (Qn , (j-')-plane (Qy be- 
ing the wave vector transfer along the chain axis) is shown 
in Fig. for r=:8.5 K. One clearly sees a "fountain" 
of magnetic scattering emanating from the Q|| = 27r/c 
ID AF zone-center. We have previously demonstrated 
that in this energy range both spin waves and continuum 
contributions to the magnetic dynamic structure factor 
are important.il However, in the present experiment sep- 
arating the two does not appear possible due to a poor 
wave vector resolution (typically 0.05 r.l.u. along the 
chains). The observed dispersion of magnetic excitations 
is consistent with the previous estimates of Vs- This is 
illustrated by the solid line in Fig. ^ that shows the calcu- 
lated Des Cloizeaux -Pearson SMA dispersion relation for 
a ID S" = 1/2 Heisenberg antiferromagnetjl^ or the lower 
bound of the 2-spinon excitation continuum,B assuming 
Vs = 130.5 meV-A. 
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FIG. 7. Neutron scattering intensity measured at 8.5 K 
using MARI as a function of wave-vector transfer along the 
chain-axis and energy transfer. The color indicates the mea- 
sured intensity according to the black and white colorbar on 
the side. The empty spaces are due to the detector arrange- 
ment which is not continuous in the 20 scattering angle. The 
black line is the lower bound of the two spinon continuum for 
J = 24.1 meV. 



D. Zone-boundary energy 

To verify that only nearest-neighbor in-chain interac- 
tions are relevant, the measured value of was compared 
to the excitation energy at the ID AF zone-boundary. 
The latter was measured at T = 1.5 K in a constant- 
Q scan at the wave vector transfer Q = (0,0,2.5) 
(Q|l = 57I-/2), using Setup IV (Fig. |). The observed 
peak was analyzed using a Gaussian profile (solid line in 
Fig. 1^) and the zone-boundary energy was determined: 
hujzB — 37.9(0.2) meV. With nearest-neighbor spin sep- 
aration along the chain axis equal to c/2, for an 5*= 1/2 
ID AF one expects htozB = 7rJ/2 and Vs = 7rcJ/4.E£l As- 
suming the nearest- neighbor model, for BaCu2Si2 07 we 
thus obtain J = 24.1(1) meV and = 130.5(5) meV-A, 
in excellent agreement with direct spin wave velocity 
measurements and magnetic susceptibility data.Q The 
agreement suggests that taking into account only nearest- 
neighbor interactions is indeed an adequate description 
of the spin chains in BaCu2Si207. 
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FIG. 8. Constant-0 scan collected at T = 1.5 K at the ID 
antiferromagnetic zone boundary Q — (0,0,2.5). The solid 
line is a Gaussian fit to the data. 



V. DISCUSSION 

Most measurements discussed in this paper have been 
performed below 4.5 rneiZ^, which is the threshold of the 
excitation continuum. eIlj Below this energy, the spec- 
trum is expected to be dominated by single-particle spin 
wave excitations. The corresponding dynamic structure 
factor, and the spin wave dispersion relation in particu- 
lar, can be calculatedjising the chain-Mean Field (chain- 
MF) approximationllj'Ej As will be discussed below, this 
theoretical construct is in remarkably good quantitative 
agreement with our experimental results on BaCu2Si207. 
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A. Model Hamiltonian 

We shall now describe the model Hamiltonian for 
BaCu2Si207 that was used as a starting point in our 
data analysis. As mentioned above, a single nearest- 
neighbor Heisenberg exchange coupling constant J is 
sufficient to describe in-chain interactions. We found 
that in order to reproduce the observed dispersion in 
the perpendicular directions, three independent inter- 
chain exchange constants are required. The Cu^"*" chains 
in BaCu2Si2 07 form a rectangular lattice with nearest- 
neighbor Cu-Cu distances along the a- and 6-axes of a/2 
and 6/2, respectively. We shall denote the corresponding 
exchange constants as Jx and Jy. In addition, we shall 
include exchange interactions along the [1,1,0] direction 
(3rd nearest-neighbor inter-chain Cu-Cu distance) in our 
model and label the corresponding coupling constant as 
J3. The model Heisenberg Hamiltonian is then written 
as: 

H = ^ ^ J Sn Sri+c/2 + JxSriSri+a/2 + JySn Sri+b/2 

i 

+ JzSriSri+a/2+b/2 + 'hS n S n+a /2-b/2- (3) 

Here the sum is taken over all spins in the system and 
is the position of spin i. 

Following Schulzll3 and treating the interchain coupling 
in a mean-field treatment, the Hamiltonian transforms 
into an effective single-chain Hamiltonian which reads 



H 



JSiS 



2NJ'ml . 



(4) 



The magnetic moment is aligned along the z-axis. Here 
N is the number of sites in the chain, mo=(— 1)' < > 
is the staggered magnetization, and h=—AJ'mo. The 
effective MF inter-chain coupling strength J' for this 
anisotropic coupling geometry is given by J' = |(2| J^^l -|- 



2 L/i, 



-4|J3|). 



In this formula S^^{Q,iij) and S'^'^{Q,uj) are dynamic 
structure factors for excitations polarized along the a- 
and 6-axes, respectively, and Q = (/i, fc, I) is the wave 
vector transfer. A is the gap induced in the quantum 
spin chains by the effective staggered exchange field in 
the magnetically ordered state. This gap corresponds to 
the excitation energy at a point in reciprocal space where 
inter-chain interactions cancel out at the RPA level, e.g., 
at Q = (0.5,0.5,1). J{Q) is the Fourier transform of 
inter-chain coupling, and is given by 



J{Q) — Jx cos(7r/i) -f Jy cos(7rfc) -|- 
J3 cos(7r(/i -I- k)) + J3 cos(7r(/i — k)). 



(8) 



Finally, the parameters Dx and Dy in Eq. ^ are 
anisotropy gaps for the two spin wave branches, that 
we empiricalljL, include in the dispersion relation, as in 
Tsukada et alu. 

The magnetic dynamic structure factor is related to 
the inelastic neutron scattering cross section through the 
polarization factors and form factors for the magnetic 
ions involved: 



dfldE' 



|/(g)|2 S--{Q,Lo)sm\Q,a)- 



SyyiQ,u)sm^iQ,b) 



(9) 
(10) 



B. Dynamic structure factor 



The model cross section defined by Eqs. p[-|l0| were 
numerically convoluted with the calculated spectrome- 
ter resolution function and used in global least-squares 
fits to the data for each experimental setting. The effec- 
tive MF inter-chain coupling strength |J'| = i(2|Jj.| -t- 
2\Jy \ + 4lJfi|) is related to the gap energy through A = 
6.175| J'|Jl3 The in-chain exchange constant was fixed at 
J = 24.1 meV, while the parameters Jx, Jy, J3, Dx 
and Dy were refined. The following values were ob- 
tained: Jx = —0.460(7) meV (ferromagnetic), Jy = 
0.200(6) meV, 2J3 = 0.152(7) meV, Dx = 0.36(2) meV 
and Dy ^ 0.21(1) meV. This yields A = 2.51 meV 
and I J'l = 0.41 meV. The obtained dispersion relation is 
shown in Fig. ^in solid lines. Symbols indicate excitation 
energies obtained in fits to individual scans. 



Esslcr et. at have previously derived the SMA dy- 
namic structure factor for weakly-coupled S = 1/2 
antiferromagnetic chains in KCuFa, using the chaip. 
MF/Random Phase Approximation (chain-MF/RPA).El 
It is straightforward to adapt their result to the coupling 
geometry for the Hamiltonian given in Eq. ^: 



5""(Q,w) oc 

syy{Q,uj) oc 



1 



LUxiQ) 

1 

UJyiQ) 



[5{u-LOx{Q))-5{uj+LOx{Q))], (5) 

[6{uj~LOy{Q))-6{Lo + LOy{Q))], (6) 



j2 sin2(7rO + 



4 

X {Jy + 2 J3 



Jy -|- 2(/3 Jx 



Jx + J{Q)) 



D^ . 



(7) 



C. Comparison with the chains-MF model 

The chain MF theory does not only reproduce the ob- 
served dispersion relation, but is in fact in good quantk 
tative agreement with the data. According to SchulzE^I 
the relations between mg, and \J'\ are given by: 



\J'\ 



Tn 



1.28A/ln(5.8J/TA,) 
mo = 1.017V|J'|/J- 



(11) 



(12) 



In our case, where T/v = 9.0 K and J = 24.1 meV, this 
gives I J'l = 0.27 meV and mo = 0.11 hb- The remain- 
ing small discrepancy with experiment can be attributed 
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to the presence of weak easy-axis (Ising-like) anisotropy, 
which clearly should have the effect of enhancing long- 
range order. Indeed, this term, no matter how small, 
induces long-range order even in the purely ID model. 

VI. CONCLUSION 

The obtained detailed characterization of BaCu2Si207 
provides a quantitative basis for the discussion of single- 
particle vs. continuum excitations in weakly-coupled spin 
chains in Zheludev et ala. The experimental data are 
remarkably well described by the quantum mean field 
model. 
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